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TeV neutrinos from microquasars in compact massive binaries 

Wlodek Bednarek 1 

ABSTRACT 

We consider a compact binary system in which a Wolf-Rayet star supplies 
matter onto a stellar mass black hole or a neutron star. This matter forms 
an accretion disk which ejects a jet as observed in Galactic microquasars. A 
part of the jet kinetic energy, typically ~ 10%, can be transfered to relativistic 
nuclei. These nuclei lose nucleons as a result of photo-disintegration process in 
collisions with thermal photons from the accretion disk and the massive star. 
Due to the head on photon-nucleus collisions most of neutrons released from 
nuclei move towards the surface of the accretion disk and/or the massive star 
producing neutrinos in collisions with the matter. We calculate the spectra of 
muon neutrinos and expected neutrino event rates in a 1 km 3 neutrino detector of 
the IceCube type from a microquasar inside our Galaxy applying, as an example, 
the parameters of the Cyg X-3 binary system. Such model predicts up to several 
neutrino events per km 3 per yr from Cyg X-3, provided that nuclei are accelerated 
to the Lorentz factors above 10 6 with a power law spectrum with an index close 
to 2. 

Subject headings: X-rays: binaries: general: individual (Cyg X-3) — stars: Wolf- 
Rayet — radiation processes: nonthermal — neutrinos — cosmic rays 



1. Introduction 

A part of X-ray binary systems, called microquasars, produce collimated outflows ob- 
served as nonthermal radio structures (e.g. Mirabel & Rodriquez 1999, Fender 2003). They 
are probably produced by relativistic jets launched from the inner parts of the accretion 
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disks around stellar mass black holes or neutron stars. High energy processes play an im- 
portant role in these sources since some of them emit high energy X-rays and 7-rays (e.g. 
Paredes 2005). However, the content of these jets, hadronic, leptonic or mixed, remains a 
puzzle. A few models have been discussed in which possible importance of hadronic pro- 
cesses in the microquasar jets can be probed by the neutrino observations. For example, 
Levinson & Waxman (2001) and Distefano et al. (2002) consider production of neutrinos in 
the interactions of protons accelerated in the jet with the radiation field of the accretion disk. 
Romero et al. (2003) discuss the case of a misaligned jet, which collide with dense wind of 
a massive star or its surface. Aharonian & Atoyan (1991) propose that relativistic hadronic 
beam can interact with blobs of matter injected by the stellar companion. Other scenarios 
for neutrino production in X-ray binaries, which do not postulate very collimated injection 
of relativistic hadrons, have been also discussed, e.g. collisions of hadrons from the pulsar 
with matter of the accretion disk (Cheng & Ruderman 1989, Anchordoqui et al. 2003) or the 
matter of the massive star (Gaisser & Stanev 1985, Kolb et al. 1985, Berezinsky et al. 1985, 
Harding & Gaisser 1990). 

In this paper we consider a different model for the neutrino production in microquasars 
in which the jet does not need to collide directly with the dense matter from the stellar 
companion (a rather rare phenomenon due to the geometrical requirements). We show that 
in very compact massive binaries containing the massive star and a compact object, e.g. 
the Cyg X-3, nuclei accelerated in the jet can lose efficiently neutrons in collisions with the 
stellar and accretion disk radiation. Most of these neutrons propagate towards the disk and 
the massive star producing neutrinos in collisions with the matter. Note that about half 
of the observed Wolf-Rayet stars, between 158 objects (van der Hucht 2001), form massive 
binary stars. Such systems should soon evolve into the binary system containing a compact 
object, i.e. similar to the Cyg X-3. In fact, a few high mass compact X-ray binaries with 
the solar mass black hole are observed at present, e.g. Cyg X-l, SS 433, V4641 Sgr (Paredes 
2005). 



2. A compact object close to a massive star 

Let's consider a solar mass black hole or a slowly rotating neutron star on a close orbit 
around a Wolf-Rayet (WR) star. The matter accreting from the massive star forms an 
accretion disk which launches a jet from its inner part as observed in microquasars. As an 
example, we consider a specific case of the Cyg X-3 binary system which basic parameters 
are the following: the Wolf-Rayet star with the surface temperature T WR = 1.36 x 10 5 K 
and radius -Rwr = l-6i2 Q , the separation of the components D = (3.6 ± 1.2)i? = 2.25/?wr- 
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Persistent X-ray emission of the accretion disk in Cyg X-3 is L d ps 10 38 erg s _1 (Bonnet- 
Bidaud & Chardin 1988). For this disk luminosity, and an assumed inner radius of the 
disk i?^ = 10 7 cm, the inner disk temperature is ~ 10 7 K. According to the the disk- 
jet symbiosis model (Falcke, Malkan & Biermann 1995), the jet power, L v is comparable 
to the observed disk luminosity, (this model has been recently successfully applied for 
different sources, see e.g. Yuan et al. 2002, Yuan, Markoff & Falcke 2002, Falcke et al. 
2004). Therefore, in farther considerations we assume L } = 0.5L d = 5 x 10 37 erg s" 1 . Note 
however, that in the outburst stages, which are relatively frequent in Cyg X-3, the X-ray 
luminosity can be several times larger. The jet in Cyg X-3 propagates with the velocity 
in the range (0.3-0.8)c (c is the velocity of light), estimated from the expension rate of 
specific components (e.g. Spencer et al. 1986) and the ratio of the jet to counter-jet flux 
(Mioduszewski et al. 2001). We apply the average value, Vj = 0.5c, reported by Marti et 
al. (2001). 

The WR stars are at the late stage of evolution of the supermassive stars which have 
already lost their hydrogen envelopes. Their winds are mainly composed from the CNO 
group of nuclei with the substantial abundance of He nuclei (see e.g. Heger & Langer 2000). 
We assume that these nuclei can be accelerated in some regions of the jet to relativistic 
energies. They create a relatively weakly relativistic blob in which nuclei with very large 
Lorentz factors are distributed isotropically in the blob reference frame. Those among these 
nuclei which happen to move in the direction of the source of soft photons, i.e. the inner part 
of the accretion disk or the disk of the massive star, suffer the most frequent disintegrations. 
As a result, most of neutrons extracted from these nuclei move along the straight lines 
towards the source of the radiation field. The number of neutrons which move towards 
the disk and the star is determined by the probability of neutron extraction from nuclei in 
such strongly anisotropic radiation defined by J A(9)fd(cos 9), where A(9) is the product of 
the cross section and the photon spectrum, / = 1 + cos 6 and 9 is the angle between the 
photon and nuclei. Let us assume that A{9) weakly depends on 9, i.e. assuming that the 
interaction regime is far away from the threshold for extraction of a nucleon from nuclei. 
Then we estimate that the probability of interaction of a nucleus which move inside the half 
hemisphere around the source of soft radiation is at least a factor 3 larger than for nucleus 
moving in the opposite half of the hemisphere. Therefore, at least 3 times more neutrons, 
extracted from nuclei, move towards the accretion disk and/or the massive star than at other 
directions. Note that at low energies the photo-disintegration cross section drops and also 
depends on /. Then significantly larger fraction of neutrons extracted from nuclei should 
move towards the source of the soft photons. 

The location of the acceleration place in the jet, the dissipation region, is unknown 
and depends on the details of the acceleration model. Therefore, we consider three different 
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possible regions inside the jet which location is determined by the strength of the radiation 
field from the disk and the massive star: 

(I) Nuclei accelerated to sufficient energies very close to the base of the jet are immedi- 
ately disintegrated in collisions with X-ray photons from the inner accretion disk. Nucleons 
from their disintegration may interact with the X-ray photons producing neutrinos through 
decay of pions, provided that they have Lorentz factors above the threshold for pion pro- 
duction, 7™ m pa / Qk-QT^ pa 10 5 , where m n is the pion mass and k-Q is the Boltzmann 
constant. Such scenario is realized if the optical depth for pion production by nucleons in 
collisions with the disk photons is above unity. The extent of the region (I) can be estimated 
by comparing the characteristic escape time of nucleons with the jet plasma, r esc = xR^/vj, 
where x is the distance along the jet in units of the inner disk radius R^, with the interaction 
time of nucleons t N£ ^ = (ndOA^c)" 1 , where a Ne = (1 — 4) x 1CT 28 cm 2 is the cross section 
for proton-photon pion production in the region of the plateau and the peak, respectively. 
The photon density from the inner disk is approximated by 

4cT SB T d 4 2 x 10 16 T 5 3 _ 3 

nd ~ q ; hp 2 ~ 2 P h - cm 1 

3c/CB-i d 2r X 2 

where T d = 10 5 T 5 is the disk temperature at its inner radius, o"sb is the Stefan-Boltzmann 
constant. Region (I) extends up to xne-^ ~ (40 — 170), i.e. at ~ (4 — 17) x 10 s cm from 
the base of the jet (see Fig. 1). We do not consider here the region (I) since it has been 
discussed recently by Levinson & Waxman (2001) and Distefano et al. (2002). 

(II) Nucleons from disintegration of nuclei are not able to produce efficiently pions in 
collisions with the disk radiation but they are able to fragment significantly on separate 
nucleons. Due to the head on collisions of the nuclei with photons, neutrons move towards 
the inner disk producing neutrinos in collisions with the matter. Such general picture (with 
the difference of injection of neutrons by protons) has been proposed for the accretion disks 
in active galactic nuclei (Nellen et al 1993). The region (II) extends along the jet to the 
distance, x^/s, at which density of photons from the disk becomes comparable to that from 
the WR star (see Fig. 1). By comparing these photon densities, we estimate x d / s pa 10 3 . 
The photon density from the massive star at the distance equal to the separation of the 
components is «wr ~ 4o"sb7wr/(3 c ^b7wr)(-Rwr/-D) 2 - For density of photons from the disk 
see Eq. 1. 

(III) In this region radiation from the WR star dominates. As we show later this region 
is at least an order of magnitude larger than the inner regions dominated by the radiation of 
the accretion disk. Neutrons extracted from nuclei in the region (III) move mainly towards 
the surface of the WR star due to the head on collisions of nuclei with the stellar photons. 



In conclusion, depending on the region in the jet where nuclei are accelerated above 
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the threshold for their efficient photo-disintegrations (I, 11, or III), we expect production 
of neutrinos either in the direction along the jet (region I), or in the direction opposite to 
the jet propagation, i.e. through the accretion disk (region II), or at the wide solid angle 
centered on the plane of the binary system (region III). Note, that in the case of significant 
Doppler boosting the borders between specific regions move down, since the disk radiation 
plays minor role due to the collimation of the relativistic nuclei in the jet. This effect is not 
considered here in detail since the Doppler factor of the jet in Cyg X-3 is not far from unity. 



3. Relativistic nuclei in the jet 

In order to estimate the maximum energies of nuclei accelerated in the jet, let us assume 
that the magnetic field energy density in the inner part of the accretion disk (at the base of 
the jet) is in the equipartition with the energy density of radiation, i.e. Ub ~ U d . For the disk 
temperature at the inner radius, T d = 10 7 K, the magnetic field can be as large as B d « 3x 10 7 
G. The magnetic field drops along the jet according to B(x) = B d /(l + ax) « B d /x (for 
x ^> I /a), where a — 0.1 rad is the jet opening angle. The acceleration rate of nuclei can 
be expressed by 

d lA /dt = £ZeBc/Am p , (2) 

where £ is the efficiency of acceleration mechanism, Am p , and Ze are the mass and the 
charge of the nuclei, and Z/A = 0.5. Nuclei can in principle reach the maximum Lorentz 
factors allowed by the spatial extend of the acceleration region, L acc ps axR d , as large as 

7^ ax « axR d (d-f A /dt)/c = £eB d R d « 10 8 £. (3) 

The maximum Lorentz factors of nuclei can be limited by their escape time along the jet, 
t csc (estimated above), or by their energy losses in collisions with thermal photons on the 
pion and the disintegration processes (the energy losses on e ± production can be neglected 
in respect to these processes). By comparing the acceleration time scale of nuclei in the 
jet, r acc = Ja/ (djA/dt) with their escape time scale, r esc , we conclude that nuclei can be 
potentially accelerated in the jet up to 7™ ax (given by Eq. 3). It has been shown above 
that in the regions II and III of the jet, the pion interaction time scale is longer than the 
escape time scale and so the acceleration time scale. Therefore, the acceleration of nuclei in 
regions II and III of the jet is not limited by their energy losses on pion production. Below 
we discuss in detail the role of disintegration process in our model. 

Heavy relativistic nuclei can lose efficiently nucleons due to the photo-disintegration 
process in collisions with thermal photons from the accretion disk or the massive star if 
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photon energies in the reference frame of the nuclei are above ~ 10 MeV, i.e. E 1 = 6/cbT7a = 
5 x 10~ 5 T 5 7a MeV > 10 MeV. This condition gives the lower limit on the Lorentz factor of 
accelerated nuclei above which their photo-disintegration becomes important, 

7 7 n « 2 x 10 5 /T 5 . (4) 

The most efficient disintegrations occur when the photon energies correspond to the maxi- 
mum in the photo-disintegration cross-section of nuclei, which is at ~ 20 MeV (with the half 
width of ~ 10 MeV) Bishop (1957), for nuclei with the mass numbers between helium and 
oxygen. 

Let's estimate the efficiency of the photo-disintegration process of nuclei in the region 
II and III. The mean free path for dissociation of a single nucleon from the nuclei is A^ 7 = 
(n<JA 7 ) _1 , where n is the density of thermal photons coming from the accretion disk (see 
Eq. 1), or the WR star, n WR (r) ks 4<t S b7wr/3c/cb^ 2 , where r = i?/i? WR , and R is the 
distance from the WR star. The photo-disintegration cross section at the peak of the giant 
resonance and the plateau region after the resonance are: a T Al xs 1.45 x 10~ 27 A cm 2 (valid 
for 7a < 3 x 10 6 ), and <r A7 ~ 1-25 x 10~ 28 A cm 2 ( 7a > 3 x 10 6 ) Karakula & Tkaczyk (1993). 
The characteristic photo-disintegration time scales of nuclei with the mass number A in the 
radiation field of the accretion disk (Td = 10 7 K) and the WR star (Twr = 1.36 x 10 5 K), 
r A7 = A A7 /c, are t£ 7 « (2.5 - 20) x 10~ 6 x 2 /A s and r^ R ~ (0.5 - A)r 2 /A s, respectively. 
On the other site, the convection escape time of nuclei along the jet is of the order of 
t c m \/R 2 — D 2 /vj, equal to ~ 8r s (if R ^> D). By comparing and r c , we get the 
upper limit for the distance from the Wolf-Rayet star, r max (2 — 15) A, below which the 
photo-disintegration process of nuclei is important. This is the upper bound on the region 
(III). Nuclei accelerated in the region (III) lose neutrons which fall onto the Wolf-Rayet star 
up at angles (3 (measured from the plane of the binary system) which fulfill the condition: 
cos/3 = D/(r max R WR ) w (0.15 - 1)/A. 

In order to determine the conditions for acceleration of nuclei let us compare the accel- 
eration time scale, r acc = 7 a/(^7a/^) ~ ^7a/(3 x 10 12 £) s, with the photo-disintegration 
time scale, ta 7 . The maximum Lorentz factor to which nuclei can be accelerated before 
significant disintegration are then: 7 ^ s ps 7.5 x 10 6 £x/A and 7 ^ = 7.5 x 10 12 ^/ (Ax), for the 
region (II) and (III) respectively (applying r = 2.25, a = 0.1, and A = 14). The comparison 
of these Lorentz factors with the minimum Lorentz factor of nuclei for the efficient photo- 
disintegration process (see Eq. 4) allows us to estimate the required value of the acceleration 
efficiency in the regions II and III, f£ w 10~ 5 (for x = 300) and f£ «3x 10~ 3 (for x = 10 4 ), 
above which nuclei at first obtain the power law spectrum and then lose separate neutrons. 
Neutrons separated from nuclei also have a power law spectrum with the spectral index of 
the parent nuclei. If £ < £ t h, neutrons are injected by nuclei with the Lorentz factors close to 
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7a"\ i- e - they are almost mono-energetic. Therefore, we consider the injection of neutrons 
by nuclei with the mono-energetic and the power law spectra. 



4. Production of neutrinos 



Most of neutrons extracted from nuclei in the region (II) move towards the inner part 
of the accretion disk and from the region (III) towards the massive star producing neutrinos 
and 7-rays in hadronic interactions with the matter. 7-rays, from the decay of neutral pions, 
are absorbed by very high column density of matter (accept a narrow region around the 
stellar limb) but neutrinos can pass through the main part of the interior of the star without 
absorption. For the parameters of the Wolf-Rayet stars the optical depth is larger than unity 
only for neutrinos moving close to the center of the star, i.e. with the impact parameter 
less than ~ 0.2i?wR which corresponds only to a few percent of the full stellar disk (see 
e.g. Bartosik et al. 2005). Therefore, the fraction of neutrinos absorbed inside the star is 
relatively low and can be neglected. We calculate the muon neutrino (and anti-neutrino) 
spectra produced by neutrons with the mono-energetic and the power law spectra. The 
multiple interactions of neutrons with the matter have been taken into account up to their 
cooling to 100 GeV. We apply the scale breaking model for hadronic interactions developed 
by Wdowczyk & Wolfendale (1987). The neutrino spectra are calculated for the mono- 
energetic nuclei accelerated in the region (II) to the Lorentz factors 2 x 10 3 and for the 
power law spectrum with the spectral indexes equal to 2 and 2.5 and the cut-offs at Lorentz 
factors 10 7 and 10 6 (see Fig. 2a). We apply the distance to the Cyg X-3 binary system equal 
to 10 kpc and assume that x — 10% of the jet energy is transported to relativistic nuclei. 
The spectra are clearly above the atmospheric neutrino background (ANB) provided that the 
spectrum of accelerated nuclei is flatter than 2.5 and extends above 7a = 10 6 . The spectra 
produced by mono-energetic nuclei are significantly below the ANB. Therefore, acceleration 
of nuclei in the region II with the efficiency £ < 10~ 5 can not be tested by neutrino telescopes. 
The neutrino spectra produced by neutrons separated from nuclei accelerated in the region 
(III) in collisions with the matter of the massive star are shown in Fig. 2b. Neutrinos 
from regions II and III have similar fluxes. They should be easily observed by the 1 km 3 
neutrino detector of the IceCube type (Hill 2001). The smaller scale detectors can already 
put interesting constraints on the efficiency of high energy processes in microquasars. For 
example, the case of mono-energetic injection of nuclei into the region (III) with the energy 
conversion efficiency x — 10% is already above the upper limit on the neutrino flux from the 
Cyg X-3 derived on the base of the 3 yr of data accumulated by the AMANDA II detector 
(Ackermann et al. 2004). This limit is also within the factor of two to the predicted neutrino 
spectrum for the case of injection of nuclei with the power law spectrum and spectral index 
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2 (see Figs. 2a and 2b). 

5. Conclusion 

Let's estimate the neutrino event rate in a 1 km 3 neutrino detector from the Cyg X-3 
type binary system. For that we apply the likelihood of neutrino detection by a 1 km 3 
detector obtained by Gaisser & Grillo (1987) and the absorption coefficients of neutrinos 
in the Earth obtained by Gandhi (2000). We consider two limiting cases, i.e. the source 
close to the horizon (neutrinos not absorbed by the Earth), and the source at the nadir 
(neutrinos partially absorbed). The number of events per year per km 3 expected for the 
mono-energetic injection of nuclei in the region III is ~ 24 events, and for nuclei with the 
power law spectrum in the region II (and the region III): 16.5 (14) events for the spectral 
index a = 2 and the cut-off at the Lorentz factor ^f^ ax = 10 7 , and 11 (8.3) events for the 
cut-off at 10 6 , and ~ 2 (~ 1) events for a = 2.5 and / y™ ax = 10 7 . Since the neutrino spectra 
do not extend much above a few 10 5 GeV, the effects of absorption of neutrinos in the Earth 
are only of the order of ~ 10%. Such event rates should be detected above the ANB by 
the IceCube detector if the spectrum of accelerated nuclei is flatter than 2.5. Note that 
predicted neutrino signals from the regions (II) and (III) should be emitted in completely 
different directions. Neutrinos from the accretion disk are collimated along the axis of the 
jet. In contrast, neutrinos from the WR star are emitted inside the broad solid angle around 
the plane of the binary system and should be modulated with the orbital period of the binary 
system reaching the maximum when the compact object is behind the Wolf-Rayet star. This 
should help to separate the neutrino signal from the ANB. Therefore, even smaller scale 
neutrino detectors, e.g. of the AMANDA II and Antares type, might be also able to detect 
such microquasars in compact massive binaries within our Galaxy or at least put strong 
constraints on the efficiency of hadron acceleration in their jets. The angular features of 
the neutrino emission, coupled with the information on the geometry of the binary system, 
should allow to distinguish between the two discussed regions for hadron acceleration giving 
information on the acceleration process and formation of jets. 

The hadronic processes considered in this paper also turn to the production of neutral 
pions which decay onto high energy gamma-rays. Gamma-rays produced in the interaction 
of neutrons with the matter of the accretion disk initiate partial cascades which electro- 
magnetic products might emerge from some parts of the accretion disk (see e.g. Cheng & 
Ruderman 1991). However, it has been shown by e.g. Carraminana (1992) and Bednarek 
(1993) that gamma-rays with energies ~ 1 — 10 3 GeV should be strongly absorbed in the 
radiation field produced by the accretion disk. Gamma-rays produced in the interaction 
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of neutrons with the matter of the massive star are absorbed inside the star except these 
ones which move inside the narrow cone around the stellar limb where the density of matter 
is low enough for the escape of gamma-rays without significant absorption in the matter. 
However these TeV gamma-rays are significantly absorbed in the thermal radiation of the 
massive star (Bednarek 2000, or Sierpowska & Bednarek 2005) since they are produced very 
close to the stellar surface. We conclude that the hadronic model considered here does not 
predict gamma-ray fluxes which can be observed by the Cherenkov telescopes at energies 
above a few tens GeV. Only gamma-rays below these energies can escape from the dense 
soft radiation of the massive star and the accretion disk being potentially detected by the 
satellite experiments such as AGILE and GLAST. However if significant amount of the jet 
energy is transfered to relativistic leptons at relatively large distances from the base of the 
jet, then microquasars can also produce observable fluxes of TeV gamma-rays in processes 
similar to that observed in BL Lacs. 

I would like to thank the anonymous referee for useful comments. This work is supported 
by the KBN grants PBZ KBN054/P03/2001 and 1P03D01028. 
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Fig. 1. — Schematic picture of a binary system containing a massive Wolf-Rayet star and a 
black hole or a neutron star (not to scale). The matter from the massive star (mainly He and 
CNO nuclei) forms an accretion disk which launches a jet. Nuclei are accelerated at the jet 
up to the Lorentz factor 7a- If their Lorentz factors are 7a > 7™"? they photo-disintegrate 
in collisions with thermal photons (e) coming from the massive star and/or the inner part 
of the accretion disk losing separate nucleons. Most of neutrons (N), separated from nuclei, 
move towards the massive star (due to the head on collisions) and/or to the inner part of 
the accretion disk and interact with the matter producing neutrinos (v). 
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Fig. 2. — The differential energy spectra of muon neutrinos, multiplied by E 2 , produced in 
hadronic interactions of neutrons with the matter of the accretion disk (a) and the massive 
Wolf Rayet star (b). The neutrino spectra are shown for the power law spectrum of nuclei, 
7^ K , described by k = 2 and the cut-off at 7™ ax = 10 7 (thick full curves) and at 10 6 (thin 
full), and k = 2.5 and 7™ ax = 10 7 (thick dashed), and for the mono-energetic nuclei with the 
Lorentz factor 2 x 10 5 (accelerated in the region III) and 2 x 10 3 (accelerated in the region II) 
(thick dot-dashed curves). The energy conversion efficiency from the jet to the relativistic 
nuclei is equal to \ = 10%. The dot-dot-dot-dashed curve indicate the atmospheric neutrino 
background within a 1° of the source (Lipari 1993), the thin dotted line shows the 3 yr 
sensitivity of the IceCube detector (Hill 2001), and the thick dotted line shows the upper 
limit obtained by AMANDA II detector for the Cyg X-3 (Ackermann et al. 2004). 



